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Abstract 

By using the quark-exchange formalism, realistic Faddeev wave functions 
and Fermi motion effect, we investigate deep inelastic electron scattering from 
A = 3 mirror nuclei in the deep- valence region. The initial valence quarks in 
put are taken from the GRV's next-to-leading order calculations on F%(x, Q 2 ) 
which give very good fit to the available data in the (x, Q 2 )-plane. It is shown 
that the free neutron to proton structure functions ratio can be extracted 
from corresponding EMC ratios for 3 He and 3 H mirror nuclei by using self- 
consistent iteration procedure and the results are in good agreement with 
other theoretical models as well as present available experimental data, espe- 
cially the one expected form the proposed 11 GeV Jefferson Laboratory. 
PACS number(s). 13.60.Hb, 21.45. +v, 14.20.Dh, 24.85. +P, 12.39.Ki 



Typeset using REVTeX 



* Corresponding author, Email : modarres@khayam.ut.ac.ir 



1 



I. INTRODUCTION 



In the framework of Standard Model, the hadrons are composed of valence- quarks, sea- 
quarks and gluons [1]. In recent years, in order to test the perturbative and non-perturbative 
nature of quantum chromodynamics, most of the experiments and theoretical works have 
been focused on studies of the hadrons (mainly proton and neutron) structure functions at 
small x ( x is Bjorken scaling variable) regions [2] where the sea-quarks and gluons play 
important roles. While the middle x region i.e. 0.3 < x < 0.7 at moderate Q 2 (the photon 
4-momentum), have been assumed to be understood and can be explained by the valence- 
quarks dynamics, there is lake of information about the nucleons structure functions in 
the deep-valence region, x > 0.7 at any Q 2 . There are a few reasons to study the quark 
distributions in nucleons at large x. The d/u quark distribution function ratio near x ~ 1 
can give information about (1) the spin-flavor symmetry breaking in the nucleons (2) the 
onset of perturbation behavior [3] and (3) search for new physics beyond Standard Model [4] 
in the high energy colliders at high Q 2 e.g. the uncertainty in the gluon distributions. 

While the proton structure functions is quite well known both experimentally and theo- 
retically [1-3], the neutron structure functions usually should be extracted from deuterium 
and because of the large nuclear corrections, there could be uncertainties as large as %50 in 
the d/u or -F 2 n /F| ~ ratios [5]. The bounds value of \ < |§r < 4 has been proposed 

u 2 

by Nachtmann [6] for the all x regions. On the other hand, based on different models the 
values of | (SU(6) symmetry) [7] , \ (phenomenological and Regge considerations) [8 ] and | 
(quarks counting roles and perturbative QCD) [9] have been predicted as x — > 1. 

Recently, the possible use of an unpolarized tritium target has been proposed, by using 
the HGeV upgraded beam of Jefferson Laboratory [10], and aimed to measure F 2 n , using 
the ratio of structure functions of helium-3 (J 7 ^ 6 ) an d tritium (J-% H ) in order to reduce 
the systematic errors both in the experimental measurements and theoretical calculations 
(which are model dependent). 

Several years ago the quark-exchange formalism was originally introduced by Hoodbhoy 
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and Jaffe [HJ] to investigate the quarks distributions in nuclear systems [11,12]. This formal- 
ism was applied by one of the authors (MM) to light nuclei [13] and nuclear matter [14] and 
was reformulated by us to derive the spin structure function of the three-nucleon systems as 
well as the proton and neutron [15] and finally it was used as the initial conditions for the 
QCD evolution equations [16] to calculate the sea-quarks and gluons contributions to the 
proton structure function at the leading and next-to-leading order (NLO) levels. In the most 
of the above calculations we found satisfactory agreements with the available experimental 
data. So we think the quark-exchange formalism is good motivation for looking at -F 2 n /.Ff 
ratios at deep- valence region i.e. x > 0.7. 

Recently, several groups have paid attention to the above matter [17-22] and have used 
different models (mostly based on impulse formalism and different spectral function approx- 
imation) to calculate F^/Ff ratios at deep- valence region. We will present their results in 
this work and compare them with ours. 

So the paper will be organized as follows: In section II we briefly explain the quark- 
exchange model and we calculate the valence quark momentum distributions for proton and 
neutron in 3 He and 3 H. In section III , we calculate the structure functions of helium 3 
and tritium by including the Fermi motion effect. The self-consistent calculation (iteration 
procedure) of the ratio / F$ will be explained in section IV. Finally in section V we present 
our numerical results, discussion and conclusion. 



II. THE QUARK-EXCHANGE FORMALISM 

Let us start with a brief summary of the quark-exchange formalism. We take the nucleon 
states to be composed of three valence quarks [11,15], 

l«) = ^ t |o) = -^AT; i ^ 3 gt ig t 2g t 3 | ) (1) 



— * — * 

where a (fXi) describe the nucleon (quark) states {P, M s , M T } ( {k, m s , m t , c})(note that 
M T {m t ) = +| and — | for the proton (up-quark) and neutron (down-quark), respectively). 
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As usual (A/" at ) denote the creation operators for the quarks (nucleons) with the state 
index /i (a). With this convention that a repeated index means a summation as well as 
integration over k. The totally antisymmetric nucleon wave function Afj? ltt2fl3 are written as, 

K m = D (vi, a*2, /is; «i) x S(h + k 2 + k 3 - P)<f>(h, h, h, P) (2) 

— * — * — * — » 

where (f>(ki, k 2 , k 3 , P) is the three nucleon wave function which is approximated by a Gaussian 
form (b ~ nucleons radius) : 

<f)(k u k 2 , k 3 ,P) = l — \ exp[-6 2 (^ 1 ^ + — ] (3) 

D(fj,i, /i 2 , A*3; ctj) are the product of four Clebsch -Gordon coefficients, C^(^ 2m , (e ClC2C3 are 
the color factors) which are defined as, 

1 1 l s l Ll s Ifl Hi 

D(fXi, V>2, ^3', OLi) — -j=€ clC2C3 —i= ^ C^ msMsa Cms^m^m 3 C^ mtMT ^ Cmf^rn^mt (4) 
Vo! V^s,t=0,l 

Now, based on the nucleon creation operators, we can define the nucleus states as, 



\Ai = 3) = (3!)-5x QlQ2Q3 AT ai V Q2 V a3t |0) (5) 

where -)/ aia2a3 are a complete antisymmetric nuclear wave functions (they are taken from 
the Faddeev calculation with Reid soft core potential [15]. According to Afnan et al. [20], 
the choice of nucleon- nucleon potential does not affect the EMC result. However we will 
investigate this matter in our future works) and it should be interpreted as the center of 
mass motion of the three nucleons. By using the same definition as the one we did for C-G 
coefficients in equation (4) i.e., 

]_ -S- --S 1.T- IIx 

D(a 1 ,a 2 ,a 3 ;A i ) = C M s l i M s M s CMl a2 M Sa3 M s C ^ 1 M T M T CMl a2 M ( 6 ) 

Then we can write the nuclear wave functions as 

X aia2a3 = X(P, Q)D( ai , a 2 , a 3 ; A,) (7) 

Relevant information comes from the momentum distribution of the constituent quarks, 
which can be defined for valence quarks with the fixed flavor and nucleon iso-spin projection 
in the three nucleon system as, 



Pmil ' Aj " (A = 3|A = 3> 



The sign bar means no summation on Mt, m t and integration over k on the repeated index 
ji. The calculation of (Ai = 3\Ai = 3) would become straightforward by doing summation 



over /x, 



(Ai = 3\Ai = 3) = ±{Ai = 3|^|A = 3) 

^,*ai«2a3 (§ ai ^ 1 § a 2^2 fiazPs _ ga 1 a2Ct 3 ,(3il3 2 (3: i ^(3 1 f32l3 3 



where 



caio^o^ AftA _ \T02 \r(3 2 \fa 3 \rf3 3 s.a 1 (3 1 /q\ 

MM M1M2M3 M2M3P1 P1P2P3 M1P2P3 V / 

After performing some algebra, one could drive the following equation for the expectation 
value of q^q, 

(Ai = 3\qtqji\Ai = 3) = 9 x «i^«3^«2« 3 A(3 2 & _ ya^c^/Ws^/Ws 

where 

w _ T2 a 3 ,/w 3 = Af£ 2a3 Af^ 2a3 6 a ^6 a ^ (10) 

and 



yaia 2 a 3 ,/3i/32/33 _ oilTai yw-a 2 yW"A yw""3 /V/7?3 

MM M°"20"3 M°"20"3 M1M2M3 P1M2M3 P1P2P3 M1M2P3 

+4A/"" 2 A/"* A/"" 3 A/"* tf Ql/31 4- 2A/"" 2 A/"* A/"" 3 A/"* I'll") 

+^- /V MMiM2- /V MM2Pi" /V PiP2P3- /V MiP2P3° + Z - /V MiM2M3- /V MM2M3- /V MM2P3- /V MiP2P3° 

Then by assuming the nucleus to be in the rest frame and defining the Fourier transform 
of x{P-i Q)i we can calculate the expectation values of different terms in equations (10) and 
(11) 



2tt 2 / " l 2 
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DQi, a 2 , a 2 , a 3 ; (h)D{a u a 2 , a 3 ; A)£>(ft, ft, ft; Ai)5 a ^5 a ^ (12) 



*a x a 2 a 3 kra x KfPi \fa 2 kT/3 2 \fa a *r/3 3 /3i/3 2 /3 3 _ 

A. ^ v ^cr 2 cT3 p.V2°3 V pip2P3 P1M2M3 P\P2P3 P\P2P3 A 

3 

/275 2 \ 5 3 - 

D(pup 2 , P3] a 3 )D(ii 1 , p 2 , p 3 ; ft)-D(«i, "2, "3; A)-D(ft, #2, ft; A) (13) 



where 



Y *«i«2a3 yyp2 yyft yy«3 yyfe xai/81 v /3i/3 2 /3 3 _ 

A J *PPlP2 J PP2Pl J P1P2P3 P\P2P3 A 

{27b 2 \^ 12 -> 
/ 1 -^j- J exp[ — -b 2 k 2 \D{p x , p 2 , p; «2)-D(pi, A*2, P5 ft)£>(Pi, P2, Ps; "3) 

£>(//!, p 2 , p 3 ; /3 3 ) J D(«i, a 2 , a 3 ; A)-D(ft, ft, ft; A)<* ai/Jl (14) 



Y «i«2«3 yw"ft yw"a3 Vfft x«i/3i /3i/3 2 /?3 

A JX Ut U9U3 J UU9U3 J U09 03 J LIT 09 03 A 



P1P-2P3 PP2P3 PP2P3 P\P2P3 

'27b r 



..b 2 \^ 

I [ J ezp[-3& A; ]D(pi, p 2 , P3; a 2 )D(p 2 , P3, p; ft)-D(p, P2, P3; "3) 



£(pi, p 2 , p 3 ; ft)^(«i, « 2 , a 3 ; A Wi, ft, ft; A)5 ai/31 (15) 



*aact 2 ct3 A/"«2 ywTfe Afas AfAs XQift v /3i/3 2 /3 3 _ 
A J P1P2P3 M2M3P1 P1P2P3 VlP2p3 A 

3 3 

Ag) ^1, ^2, A*3| a 2 )D(p u p 2 , p 3 ; ft)D(pi, p 2 , p 3 ; a 3 ) 
L>(p!, p 2 , p 3 ; ft) J D(a 1 , a 2 , a 3 ; A)£>(ft, ft, ft; A)5 aift (16) 



/•oo ;-oo ^1 3X 2 

/ = 8?t 2 / x 2 dx y 2 dy / d(cos#)e:rp[ — — ]|x(a;, cos#)| 2 (17) 
Jo Jo J-i 4tr 

The above equations have been calculated with the same approximation as the one used in 

the references [11,16], specially a leading order expansion for x(PiQ) [16]. AS we pointed 

out before according to references [20,23], the other choice of nucleon-nucleon potentials do 

not change the ^4 = 3 nuclear wave function and the EMC effect very much . 
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III. NUCLEUS STRUCTURE FUNCTION 



The structure function measures the distribution of quarks as a function of k + (the light- 
cone momentum of initial quark) in the target rest fame which is equivalent to boosting the 
nucleus to an infinite momentum fame. This is usually done by using an ad hoc prescription 
for k° as a function of \k\(k° = [(k 2 + m 2 )s — e ]). It has been shown that the resulting 
structure functions are not sensitive to this assumption [11,24]. So the valence-quark distri- 
bution at each Q 2 , can be related to momentum distribution for each flavor in the nucleons 
of nucleus Ai according to the following equation (j — p,n (a — u, d) for protons (up-quarks) 
and neutrons (down-quarks), respectively), 

qj a (x, Q 2 ; A) = J pi(k; A)5{x - ^±) dk (18) 

After doing the angular integration, we get, 

poo _, 

qj a (x, Q 2 ; Ai) = 2ttM / pi(k; Ai)k dk (19) 

with 

k (x) - {XM + eo)2 " ^ (20) 

where m (M) is the quark (nucleon) mass and eo is the quark binding energy. For each Q 2 
value, it is possible (by using the fitting procedure as it will be described later on) to calculate 
the corresponding values of m and eo- Finally, the target structure function F^(x, Q 2 ) can 
be expressed in terms of the valence quark distributions as following : 

F 2 %(x,Q 2 )=x J2 Q 2 aq] a (x,Q 2 ;A t ) (21) 

a=u,d;j=p,n 

So we have found the relation between quark momentum distribution and the target struc- 
ture function. In order to fix the values of m and eo, we apply the above equations as well as 
the quark-exchange formalism to the proton (Ai = I, Mr = M-t = \) as our target (obvi- 
ously by considering the proton as our target there is no exchange term and we have just the 
direct term). In this case for each value of b we find the pairs m and eo such that we get the 



best fit to the valence u and d quarks distribution functions of GRV's partons structure func- 
tion calculations [25]. The GRV's partons structure functions fit the experimental proton 
structure function data over the whole range of (x, Q 2 ) plane very well. Figure 1 shows our 
fitted proton structure function (obviously only for valence quarks) for b = 0.8/m with the 
(m,e ) pairs of (165MeV, 150MeV) and (165MeV, 250MeV) at Q 2 = 4GeV 2 . The dotted 
curve is those of GRV and other curves have been obtained by the same pairs of (m, eo), as 
above, but for different b. The experimental data of SLAC [19] and NMC [29] as well as 
the full GRV's NLO structure functions for proton (dash curve) and neutron (dash dotted 
curve) have been also given for comparison. It is seen that for x > 0.2 we get very good fit 
to the data. With the same parameters we can calculate the corresponding neutron, 3 He 
and 3 if structure functions which are uncertain because of the lake of information about the 
neutron structure functions. Obviously in general our method fails for x — > (because of 
fitting procedure). It is also not good as x — > 1 especially for the nucleus target, because we 
have ignored the Fermi motion effect by using the expansion in equation (17). In order to 
take into the account the Fermi motion effect we have worked in the convolutions approach 
and the harmonic oscillator basis for tritium and helium 3 with the procedure described in 
reference [1,26] as following. In the convolutions approach the nucleus structure function 
can be written as : 

F^(x,Q 2 )= W dzf^(z)Fi(x/z,Q 2 ) (22) 

j=p,n Jx 

where /^(z) and F£(x,Q 2 ) are the proton or neutron distribution functions in the target 
nucleus and the corresponding structure functions, respectively. In the harmonic oscillator 
basis we have: 

//*(*) = E Gn J ,l j <S nj ,l j (z, M, kuj, Srijlj) (23) 
rij,lj 

ft 2 2 

where G n j, S n j(z, M,hu,e n i), hu = and e n i are the occupation numbers, sum of har- 
monic oscillators polynomials (for present calculation it is just a Gaussian function), oscil- 
lator parameter and single particle energies [26]. For three-body system a 2 = 2< |j a> where 
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the rms radius ,(< r 2 >)5, is 1.95 fm and 1.7 /m for helium 3 and tritium [27], respectively. 
We also set e n i = 0, since we only intend to calculate the Fermi motion effect (the quark 
exchange is responsible for the binding effect). Figure 2(a) shows the ratio of Fermi motion 
effect for 3 He to 3 H i.e TZ 3 F H J mi {x, Q 2 )jlt^ erm ^x, Q 2 ) where, 

q-, 3 He, Fermi/ _ Fermi,lk X i Q ) sp 3 H, Fermi/ _ ^ ~ Fermi^k X i Q ) 

EMC [X > * j " 2F$(x, Q 2 ) + F?(x, Q 2 ) ' EMC } ~ F£(x, Q 2 ) + 2F 2 "(x, Q 2 ) ' 

(24) 

It is seen from this figure that the Fermi motion is the same for helium 3 and tritium up to 
x ~ 0.55 and since the size of 3 He is larger than 3 H, the ratio start to decrease from this 
value i.e. the Fermi motion effect has larger size in tritium with respect to helium, as one 
expects. 

Finally the total structure function for each nucleus can be written as the sum of the 
quark-exchange, equation (21) and Fermi motion, equation (22), effects : 

Ff*(x, Q 2 ) = F 2 %(x, Q 2 ) + F 2 % rmi (x, Q 2 ) (25) 

We hope in our future works we omit this approximation by calculating the Fermi motion 
effect in the framework of quark-exchange formalism and the Faddeev wave functions for 
three-nucleon systems. 

IV. SELF-CONSISTENT TREATMENT OF F?/Ff 

By defining the EMC-type ratios for the structure functions of helium 3 and tritium (as 
the one we did for the Fermi motion effect) i.e.: 

Kmc(x, Q 2 ) = 2F?faQ*)+%lx,Q*) (26) 

and 

* c(x ' Q2) = F 2 p (x,Q 2 )[2%\x,Q 2 ) (27) 
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we can calculate the above EMC ratios by using the fitted GRV's proton and neutron 
structure functions, the convolutions approach and the quark-exchange formalism. In figures 
2(b) and 2(c) we have plotted these ratios for different values of b. The HERMES helium 3 
data is taken form reference [28] at Q 2 ~ 7GeV 2 . We should point out here that the quoted 
data is the combination of helium 3, deuterium and proton cross-sections (see Ackerstaff 



[28]) (7Z = r d +F p ) ■ We get good agreement with the present available data, especially for 
the deep-valence region. It is interesting that in this region the EMC ratios are not very 
sensitive to the different values of b. 

Now by dividing the above two EMC ratios we find the following equation: 

K He/3H (x O 2 ) - n *™cM 2 ) _ n >He/>H (x Q 2 l + 2F?(x,Q 2 )/F?( x ,Q 2 ) 

with 

K -»«, v , Q > ) = ^g (29) 

and by imposing the constrain, 

n E H M'c H M 2 )^^ (so) 

which will be discussed later on. 

The above equation can be solved for the neutron to proton structure functions ratio, 
F 2 n (x,Q 2 )/F 2 p (x,Q 2 ), in terms of the EMC ratio, U^ 11 (x,Q 2 ), which directly yields : 

F?(x,Q 2 ) = 2n E H ^ H (x,Q 2 )-K iH */ :iH (x,Q 2 ) 
F%(x,Q 2 ) 2n 3 ^(x,Q 2 )-n i E H ^ H (x,Q 2 ) 

Equations (28) and (31) are coupled in terms of the neutron to proton structure function 
ratio and will be solved by a self-consistent iteration procedure i.e the out put of equation 
(31) will be used as the in put for equation (28) and so on, by using equation (30) as our 
constrain. 
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V. RESULT, DISCUSSION AND CONCLUSION 



In figure 3 the ratios of structure functions of helium 3 to tritium (equation (29)) have 
been plotted for different values of b (full curves) by using GRV's structure functions. The 
sun-burst points are the expected ratio which have been estimated by using the kinematics 
of the proposed 11 GeV Jefferson Laboratory experiment [10,20]. It is seen that there is very 
good agreement between the estimated prediction [10,20] and the present calculation. This 
shows that the GRV's proton and neutron structure functions and present model without 
any parameter can predict reasonabley the structure functions of helium 3 and tritium. 

The EMC ratios (equation (28)) has been given in figure 4. In general, since we have 
taken into the account all of the properties of the structure functions of proton, neutron, 
helium 3 and tritium, it is expected that TZ^ C H {x, Q 2 ) ~ 1, as we imposed it as a constrain 
in equation (30). However it is seen that the calculated EMC ratio has a very small variation 
from 1Zem C H (x,Q 2 ) = 1. But this deviation increases as x — > 1. The similar behavior is 
also seen in the others calculations in which the impulse approximation have been used [20]. 

Finally, in figure 5 the calculated neutron to proton structure functions ratios have 
been plotted for different values of b. The curve without iteration is the one in which we 
have not imposed the constrain equation (30). The first and second iteration curves have 
been calculated by considering the constrain TZ^q H (x, Q 2 ) = 1.01. The third iteration 
is not distinguishable from the second one. The limit values of this ratio for x — 1 (as it 
was discussed in the introduction) from different models are shown in the left part of this 
figure by arrows. The data are from references [17,18,19,28]. The iterations approximately 
converge after the third iteration. Our results cover all of the range of present available 
data. This can be done by varying TZ^ C H {x,Q 2 ) between 0.99 and 1.01. Variation of b, 
i.e. the nucleon size, does not affect the results. But the size of EMC ratios is important 
as x becomes closer to 1. Our calculated neutron to proton structure functions ratios are 
also in good agreement with present theoretical calculation in which different models and 
approximations have been used [17-22]. 
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In conclusion, we have calculated the neutron to proton structure functions ratios in the 
frame-work of quark-exchange model. We have treated u and d quarks as well as proton and 
neutron explicitly in our formalism and we have found satisfactory results compare to present 
available data and others theoretical models. However, we can improve our calculation by (i) 
treating the Fermi motion effect explicitly in the quark-exchange frame-work i.e. calculating 
the expansion we made in equations (13)-(16) as well as using other choice of nuclear wave- 
functions which have been calculated with the new nucleon-nucleon potentials, (ii) evaluating 
the connected three-body diagram which we have been ignored in the present calculation 
[15] and finally (iii) taking into the account the nucleon-nucleon correlations. The later 
is interesting, because for nuclei one could have non-vanishing structure functions for x 
larger than one [29,30], which can be measured and calculated by considering the effect of 
short-range correlations in nuclei. 

MM would like to thank the University of Tehran for supporting him under the grants 
provided by its Research Council. 
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FIGURES 

FIG. 1. The proton structure function, only with valence quarks, for (m, eo) pairs of 
(!Q5MeV, 150MeV) and (165MeV, 250MeV) at Q 2 = AGeV 2 and various b values. The pairs 
(165MeF, 150MeV) and (165MeF, 250MeV) at Q 2 = AGeV 2 have been chosen such that 
b = 0.8 fm gives the best fit to the GRV's valence-quarks distributions (dotted curve) [25]. The 
dash and dotted dash curves are the GRV's full NLO proton and neutron distribution functions. 
The NMC and SLAC data [19,28] for the structure function of proton have been also given for 
comparison. 

FIG. 2. (a) The EMC ratios of the structure functions of helium 3 to tritium by considering 
only the Fermi motion effect with parameters given in the text. The EMC effect in helium 3 (b) 
and tritium (c) for different values of b. 

FIG. 3. (a) The ratios of structure functions of helium 3 to tritium (equation (29))for different 
values of b (full curves). The dash curve is after second iteration with TZ^^ C H (x,Q 2 ) = 1.01 . 
The sun-burst point are from references [10,20] (see also the text). 

FIG. 4. The EMC ratios for different values of b without any constrain (equation (28)). 

FIG. 5. The neutron to proton structure functions ratios (equation (31), see the text for more 
details ). The data are from Whitlow et al. [19], Melnitchouk and Thomas [18], Bodek et al. [17] 
and NMC [28]. 
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